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Abstract 

The pion source as seen through HBT correlations at RHIC energies is investigated within the 
UrQMD approach. We find that the calculated transverse momentum, centrality, and system size 
dependence of the Pratt-HBT radii Rl and Rs are reasonably well in line with experimental data. 
The predicted Ro values in central heavy ion collisions are larger as compared to experimental 
data. The corresponding quantity \J Rq — R 2 S of the pion emission source is somewhat larger than 
experimental estimates. 
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In the quest to discover the high temperature phase of Quantum Chromodynamics 
(QCD), the Quark Gluon Plasma (QGP), the beam energies of accelerators have been 
boosted upwards from SIS, AGS, SPS, to RHIC. However, it is well known that the phase 
transition from hadrons to quarks might only occur in a small volume part of the system 
and within a rather short timespan in heavy ion collisions (HICs). This implies that the 
QGP drops formed at RHIC might be represented only by a few, locally thermally equili- 
brated drops of matter, in which quarks and gluons are de-confined. Thus, it is essential to 
probe the space-time structure of the (equilibrated?) source - the "region of homogeneity". 
Unfortunately, the small size and transient nature of the reactions preclude direct measure- 
ment of the time and/or position. Instead, correlations of two final-state particles at small 
relative momenta provide the most direct link to the size and lifetime of subatomic reac- 
tions. A well-established technique, the so-called "femtoscopy" or "HBT" in the heavy- ion 



community (in reference to Hanbury-Brown and Twiss's original work withphotons) 
been extensively used for HICs with energies from SIS, AGS, SPS, to Rt 
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Numerous HBT-measurements with various two-particle species have been pursued (see 



e.g. 



If] | and references therein). Identically charged pion interferometry has been most 



extensively investigated. Basic, but important systematics o 



femtoscopic measurements 



from the AGS, SPS, and RHIC have been discovered [Jj, 13, [1^, [16 ] , such as the dependence 
of the HBT radii on system size, collision centrality, rapidity, transverse momentum, and 
particle mass. However, the existence of the so-called HBT-puzzle (i.e., the fact that model 



calculations that incorporate a phase transition to a new state of matter with 



of freedom significantly over-predict the observed source sizes) jl|, |4 , 
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drives us to a deeper and more systematical theoretical exploration. The Ultra-relativistic 



Quantum Molecular Dynamics (UrQMD, v2.2) transport model (employing hadronic and 

|) and the 



string degrees of freedom) (for details, _the_ reader is referred to Refs. 
analyzing program CRAB (v3.0/5) 
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particle interferometry. It is known that the probable QGP phase is not strictly treated in 
the present version of UrQMD model for the early-stage HICs. However, the failure of fluid 
dynamical models to get the /c^-dependence of the HBT radii jjl] suggests that flow is not the 
only aspect that influences the observed /^-dependence of HBT radii. A realistic hadronic 
cascade model such as UrQMD can thus throw light on what other mechanisms during the 
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late freeze-out stage induce strong coordinate-momentum correlations and thereby generate 
the observed strong fcy-dependence of the HBT radii. With this equipment, the excitation 
functions of the HBT radii of negatively charged pions are calculated systematically. In 
this paper, we focus on RHIC energies, where the biggest challenge is faced by the current 
theoretical models. The femtoscopy results at lower energies will be presented in a further 



study 



4l|. 



The correlation function of two particles is decomposed in Pratt's (so-called longitudinal 
co-moving system "Out-Side-Long") three-dimensional convention (Pratt-radii). The three- 
dimensional correlation function is fit with the standard Gaussian form: 



C(qo, Qs, Ql) = 1 + Xexp(-R 2 L q 2 L - R 2 q 2 - R 2 s q 2 s - 2R 2 OL q q L ), (1) 

in which qi and Ri are the components of the pair momentum difference q = P2 — Pi 
and the homogeneity length (Pratt-radii) in the i direction, respectively. The A is the 
incoherence factor, which lies between (complete coherence) and 1 (complete incoherence) 
in realistic HICs. Rol represents the cross-term. For mass-symmetric colliding system, 
the Rol vanishes automatically at mid-rapidity due to the longitudinal reflection symmetry 
and is also found to be negligible in our present calculations. Furthermore, in the present 
UrQMD calculations at RHIC energies, the Coulomb and other potential interactions are 
not considered (the "cascade mode" is used) due to the excessive computing times which 
would have been used otherwise. The Coulomb final state interactions are not taken into 
account in the analyzing program CRAB. 

Fig. Ogives the transverse momentum dependence (k^ = (p 1T + p 2T )/2) of the Pratt- 
radii Rl (left plots), Ro (middle plots), and Rs (right plots) at nucleon-nucleon center-of- 
mass energies ^/snn = 30, 62.4, 130, and 200 GeV (plots from top to bottom) in Au+Au 
reactions. The experimental results at energies -y/sjviv = 62.4, 130, and 200 GeV for central 
collisions (< 15%, < 10%, and < 5% of the total cross section <7t, respectively) and at 



mid-rapidities (\rj C m\ < 0.5) are taken from Refs. |22l. l23l. l24l. l25l. l26l| . The experimental error 
bars are shown as the sum of both statistical and systematic errors. The corresponding 
calculations with the same trigger- and acceptance- conditions as in the experiments are 
shown, as well as the lower energy case ^/s^n = 30 GeV for central collisions (< 15% of 
<t t ). 

Both the absolute values and the decrease of the Pratt-radii Rl and Rs with transverse 



momentum is reproduced by the present model calculations very well. The origin of the 
decrease of the Pratt-radii with the increase of transverse momentum is still under discussion: 
it may be caused by the strong underlying transverse how Q, or, by the temperature 
inhomogeneities within the hadron source (point of view of the hydrodynamics model) [31] . 
Here, it is also seen that the calculated fey-dependence of Rs is somewhat flatter than that 
of Rl, which implies that flow effects on the fey- dependence of the Pratt-radii can at least 
not be excluded. Besides the flow effect, the surface-like emission charactistic of microscopic 
models should play significantly (or even dominant) role on HBT parameters as well because 
also other Cascade/Boltzmann calculations (see e.g., the Relativistic Quantum Molecular 
Dynamics model (RQMD) , the Hadronic Rescattering Model (HRM) |3j, and A 

Multi-Phase Transport model (AMPT) 4||) can reproduce the fey dependence of Pratt radii 
(almost) equally well. The transport of the hadronic gas, or the final state hadronic 'corona' 
dynamics, should be further investigated since the QGP image is obscured by the hadronic 
'corona' ^3, 42, 3,^^. The UrQMD calculations of Rl and Rs reproduce the experimental 
data well within the error bars, while the calculated i?o's are larger than the experimental 
data — the Ro is about 25% too large. We must conclude that at RHIC, larger ratios of Ro 
and Rs are seen from hadron transport model than expected. Similar observations have also 
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been reported from most of other model calculations (c.f. 
Ref. |28J] has argued that the origin of this HBT-puzzle might be multifaceted. 

Fig. 121 shows the fey-dependence of the Pratt-radii in Au+Au reaction at ^/snn = 200 
GeV for four centralities: - 5%, 10 - 20%, 30 - 50%, and 50 - 80% of total cross section. 
Here, a pseudo-rapidity cut \r] cm \ < 0.5 has been chosen. For better visibility, we have 
shifted in the figure the values of the radii by 0, 5, 10, and 15 fm for the four centralities. It 
is very interesting to see that our calculations for the centrality dependence of Pratt-radii 
are in reasonably good agreement with the experimental data. The important observation, 
however, is that the calculated Ro values tend to deviate from the data for central reactions 
by ~ 20%, while they agree at midcentral and peripheral collisions. 

The centrality dependence of the Pratt-radii can be seen more clearly from Fig. 01 , which 
shows the Pratt-radii at fey = 250 ~ 350MeV/c as a function of the number of participants 



Npart- The quantity W R Q — R s is also shown for comparison. In spite of the reasonable 



results on the centrality dependence of the Pratt-radii, the calculated quantity y R Q — R 2 S 
obviously deviates from that extracted from data for the most central collisions: it is about 
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FIG. 1: (Color online) Transverse momentum kx dependence (at midrapidity) of the Pratt-radii 
Rl, Rqi and -R^in Au+Au collisions at -v/sat /v = 30, 62.4, 130, and 200 GeV. Experimental data 
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261 ] . The experimental errors are the sums 



for the latter three cases are also shown 
of both statistical and systematic errors. 

twice as large as measured by experiments. In previous UrQMD calculations on elliptic flow 
at RHIC it was found that only ~ 60% of the observed elliptic flow is produced [37[, which 
is probably due to a smaller anisotropy in the pressure gradients in the early stage of RHIC 
collisions in the UrQMD simulations compared to hydrodynamics. Based on Ref. [13], the 
Ro contains the contributions from temporal extent of the source and becomes larger with a 
smaller transverse freeze-out momentum of particle pairs with a certain duration time. While 
the expansion has no effect on the Rs- The calculated larger R Q and, correspondingly, the 



larger calculated quantity y R<j — R% m comparison to the data might therefore be related to 
the elliptic flow problem in the model. The pion freeze-out volume Vf has been investigated 
thoroughly experimentally j2|. Vf can be expressed as Vf = (27r) 3//2 i?ii?|. The linear 
increase of Vf with N part is expected if the pions freeze out at a constant density pf at a 
certain beam energy as observed in Ref. [?J and implied in Ref. jlj]. It can be explained 
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FIG. 2: (Color online) Midrapidity /^-dependence of Pratt-radii in Au+Au reaction at ^snn = 200 
GeV for four centralizes: - 5%, 10 - 20%, 30 - 50%, and 50 - 80% of the total cross section, 
which are shifted by 15, 10, 5, and fm, respectively. Experimental data are taken from Ref. 124]. 



reasonably well by the present model, although a smaller "thermal ellipse" is predicted due 
to a little shorter and Rs values, as shown. 

The calculations for central Cu+Cu collisions are shown in Fig. |3] (solid dots), which 
are in line with the centrality dependence of the HBT space-time structure calculated for 
Au+Au collisions. This implies that the participant multiplicity is a very good scaling 
variable, which drives the geometry (HBT radii) at mid-rapidity, at least for mid-size to 
heavy systems. In order to check this, the fey-dependence of the ratios of the Pratt-radii 
between different systems is shown in Fig. 0]. The radius ratios shown are from (a), Cu+Cu 
vs. p+p, (b), Au+Au vs. p+p, and (c), Au+Au vs. Cu+Cu. In order to read the 
figure more conveniently, the Ro and Rs ratios are shifted by 5 and 10, respectively. In 
the p+p calculation, the non-femtoscopic correlations at large relative momenta are also 
seen, that is, the pion correlations function saturates at large relative momentum, but the 
value is not equal to 1, which was also implied by the preliminary data reported recently 
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FIG. 3: (Color online) Centrality dependence of the Pratt-radii ((a)-(c)), the quantity J Rq — R s 
(in (d)), and the freeze-out volume Vf ((e)) at hr = 250 ~ 350MeV/c, in Au+Au collisions at 
yj $nn = 200 GeV. Experimental data are taken from Ref. 2^|. The calculated results for central 
Cu+Cu collisions shown with solid dots are perfectly located on the Au+Au systematic curves. 

by Ref. Q. We eliminate this effect by multiplying a constant into the parametrization 
of the correlation function. The Rl and Rs values in p+p collisions can be reproduced 
well, while the calculated Ro values are again larger than the experimental data, similar to 
the nucleus-nucleus collisions. This might be the origin of the whole puzzle, namely that 
the HBT-correlations are somewhat incorrectly put into the model in the elementary p+p 
dynamics. Figs. |U (a) and (b) show that the calculated Rl and Rs ratios, reproduce the 
experimental data reasonably well. They are almost flat as a function of k?, which means 
the ^-dependence of the Pratt-radii still exists in the elementary p+p collisions at RHIC 
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FIG. 4: (Color online) /cr-dependence of the ratios of Pratt-radii between different systems: (a), 
Cu+Cu vs. p+p, (b), Au+Au vs. p+p, and (c), Au+Au vs. Cu+Cu. Central collisions at 
\J $nn = 200 GeV and a mid-rapidity cut are chosen. In (a) and (b), the ratios are shifted by 5 
and 10, respectively. The preliminary experimental data are taken from Ref. jl^J ]. 

energies. 

The fact that the kx dependence of Rl and Rs radii in p+p collisions is similar to AA 
reaction is puzzling at first glance. In order to explore the origin of the kx dependence in pp 
and AA, we randomly exchange the momentum vectors of the pions at freeze-out and recal- 
culated the fey-dependence of Pratt radii in the p+p (solid lines with open square symbols) 
and the most central Au+Au (dotted lines with open square symbols, rescaled by the Pratt 
radii at kx = 200 MeV/c in the p+p collisions) collisions at RHIC energy ^snn = 200 
GeV, as shown in Fig. El We also show the standard calculation results (solid squares) 
and the experimental data (stars) in p+p collisions. After considering the random mixture 
of the momenta of freeze-out pions, the kx dependence of Pratt radii essentially vanishes, 
especially in the transverse direction. This is a clear indication that space-momentum corre- 
lations drive the behavior of the Pratt radii with increasing k? both in AA and pp! However, 
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FIG. 5: (Color online) A?r-dependence of the Pratt-radii Rl (left plot), Ro (middle), and Rs (right) 
in the p+p and the most central Au+Au collisions at -^/saT/v = 200 GeV. The Pratt radii in Au+Au 
collisions are rescaled by the Pratt radii at kj- = 200 MeV/c in the p+p collisions. The standard 
calculations (titled as " standard" ) as well as the calculations after considering the random mixture 
of the momenta of freeze-out pions (titled as "p-mixed") are shown. The preliminary experimental 
data for p+p collisions are shown with stars. 

it is important to stress that the origin of the space-momentum correlation in p+p is most 
probably due to jet-like structures and not flow. 

Since the UrQMD model gives too large a Ro value in p+p collisions, the calculated 
Ro-r&tio between Au+Au (or Cu+Cu) and p+p is smaller than the experimental data, in 
particular at low transverse momenta. This phenomenon disappears when we consider the 
Pratt-radii-ratios between two heavy systems, for examples, between Au+Au and Cu+Cu, 
in Fig. El (c). It is interesting to see that for the Pratt-radii-ratios between Au+Au and 
Cu+Cu collisions, all radii-ratios are flat with k T and approach ~ 1.4, which is equal to the 
ratio between the initial radii of nuclei. 

To summarize, by using the CRAB program, we analyzed the evolution of the Pratt- 
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radii Rl, Ro, and Rs at RHIC energies in collisions simulated by the UrQMD transport 
model. The calculated transverse momentum-, centrality-, and system dependence of the 
Pratt-radii are in reasonable agreement with the experimental data. The calculated Ro 
values for central collisions are ~ 25% larger as compared to experimental data. As a 
consequence, the extracted quantity yj Rq — i?| of the pion emission source is somewhat 
larger than experimental estimates. 
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